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ABSTRACT 

We introduce a metallicity dependence of Type la supernova (SN la) rate into the Galactic and cosmic 
chemical evolution models. In our SN la progenitor scenario, the accreting white dwarf (WD) blows a 
strong wind to reach the Chandrasekhar mass limit. If the iron abundance of the progenitors is as low 
as [Fe/H] < — 1, then the wind is too weak for SNe la to occur. Our model successfully reproduces 
the observed chemical evolution in the solar neighborhood. We make the following predictions which 
can test this metallicity effect: 1) SNe la are not found in the low-iron abundance environments such 
as dwarf galaxies and the outskirts of spirals. 2) The cosmic SN la rate drops at z ^ 1 — 2 due to the 
low-iron abundance, which can be observed with the Next Generation Space Telescope. At z > 1 — 2, 
SNe la can be found only in the environments where the timescale of metal enrichment is sufficiently 
short as in starburst galaxies and ellipticals. 

The low-metallicity inhibition of SNe la can shed new light on the following issues: 1) The limited 
metallicity range of the SN la progenitors would imply that "evolution effects" are relatively small for 
the use of high redshift SNe la to determine the cosmological parameters. 2) WDs of halo populations 
are poor producers of SNe la, so that the WD contribution to the halo mass is not constrained from the 
iron abundance in the halo. 3) The abundance patterns of globular clusters and field stars in the Galactic 
halo lack of SN la signatures in spite of their age difference of several Gyrs, which can be explained by 
the low-metallicity inhibition of SNe la. 4) It could also explain why the SN la contamination is not 
seen in the damped Lya systems for over a wide range of redshift. 

Subject headings: Abundances — binaries: close — Cosmology: general — Galaxy: evolution — 
galaxies: evolution — stars: supernovae 



1. INTRODUCTION 

There exist two distinct types of supernova explosions: 
One is Type II supernovae (SNe II), which are the core 
collapse-induced explosions of short-lived massive stars 
(> 8Mq) and produce more O and Mg relative to Fe (i.e., 
[0/Fe] > 0), and the other is Type la supernovae (SNe 



tain progenitor systems of SNe la; for example, Yoshii, 
Tsujimoto, & Nomoto (1996) estimated the lifetime of SN 
la progenitors as long as 0.5 — 3 Gyr. 

SNe la have been discove red up to z ^ 0.93 by Su- 



pernova Cosmology Project ( Perlmutter et al. 1997) an d 
High-z Supernova Search Team (Garnavich et al. 1998) 



la), which are the ihcrnionuclear explosions of accreting 
white dwaifs (WDs) m close buiaiies and piuduce uiusLl^ 
Fe and little O. The exact companion stars of the WDs 
have not been identified but must be relatively long-lived 
low mass stars. 

The role of these two types of supernovae in the chemical 



They hav e given the SN la rate a.t z ^ 0.4 (Pain et 



al. 1996 ) but will provide the SN la rate history over 
< z < 1. With the Next Generation Space Telescope, 
both SNe la and II will be observed through z 4. In 
theoretical approach, the cosmic SN la rate as a function 
of redshift has been constructed for a cosmic star forma- 



evolution of galaxies can be seen in the abundance ratios 



uf ytfiiy wlLh dllltiiuiiL luuLstUluUluy, mum iiuLubli; lii Lliu 

[0/Fe]-[Fu/II] lelaLiuu. Metal puui sLais with [Fc/II] < -1 



tion rate (SFR) ( Ruiz-Lapuente fc Canal 199S Yungelsoi 
fc Liyio 1998| ; [Sadat et al. 1998| ; |Madau, Valle fc Panagia 



1998). The comparison between the model prediction and 



have [0/Fe] ^ 0.45 on the average (Nissen et al. 1994 



Gratton 1991), while disk stars with [I'e/HJ > — 1 sho' 



a decrease in fO/Fe] with increasing metallicity 
' 3|; iBarbuy fc Erdelyi-Mendes 198"9| 



son et al. 199S 




observations can constrain the lifetime of the progenitor 
systems of SNe la. 

The progenitors of the majority of SNe la are most 
likely the Chandrasekhar (C h) mass WDs (e.g., INomoto 



Iwamoto fc Kishimoto 1997a for a recent review) , although 



ton 1991). Such an evolutionary change in [0/FeJ against 
[Fe/HJ has been explained with the early heavy element 
product ion by SNe II and the delay ed enrichment of Fe by 
SNe la (Matteucci & Greggio 1986). Conversely, chemical 
evolution models can constrain the nature of still uncer- 



the sub-Ch mass models might correspond to some pecu- 
liar sub-luminous SNe la. The early time spectra of the 
majority of SNe la are in excellent agreement with the syn- 
thetic spectra of the Ch mass models, while the spectra of 
the sub-Ch mass models are too blue to be compatible 
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with obse rvations ( Hoflich fc Khokhlov 1996 ; Nugent 
al. 1397). For the evolution of accreting WDs toward 



the Ch mass, two scenarios have been proposed: One is 
a double degenerate (DD) scenario, i.e., merging of dou- 
ble C+O WDs with a combined mass surpassin g the Ch 
mass limit (jlben fc Tutukov 1984|; IWebbink 19841), and the 



other is a smglc degenerate (SD) scenario, i.e., accretion 



of hydrogen-rich matter via mass tra nsfer from a binary 
companion (e.g., Nomoto et al. 1994| for a review). The 



issue of DD vs. SD is still debated (e.g., Branch et al. 1995 



for a review), although theoretical modeling has indicated 
that the merging of W Ds does not make typical SNe la 
dSaio fc Nomoto 1985| , 1998). 

For the SD scenario, a new evolutionary model has been 
proposed by Hachisu, Kato, & Nomoto (1996, 1998; here- 
after IIKN96 and HKN98, respectively) and Hachisu & 
Kato (1998; hereafter HK98). HKN96 have shown that if 
the accretion rate exceeds a certain limit, the WD blows 
a strong wind and burns hydrogen steadily to increase the 
WD mass (see section 2). HKN98 have further invoked 
the effect of stripping-off of the envelope from the com- 
panion by the strong wind, and shown that the WD can 
reach the Ch mass for much wider binary parameter space 
than found by HKN96, Li & van den Heuvel (1997), and 
Yungelson & Livio (1998); the allowed parameter space 
may be large enough to account for the SN la frequency. 
Moreover, HK98 have found an important metallicity ef- 
fect; if the iron abundance of the accreted matter is as low 
as [Fe/H] < — 1, the WD wind is too weak to increase the 
WD mass through the Ch mass. 

In the present Letter, we apply the above two scenarios 
to the chemical evolution models, and compare the cases 
with and without the metallicity effect on SNe la. We have 
found that the model for SD scenario with the metallicity 
effect is significantly better to reproduce the evolutionary 
change in [0/Fe] and other properties (Section 3). Using 
the metallicity dependent SN la rate, we make a prediction 
for the cosmic supernova rate history (Section 4). In sec- 
tion 5 we discuss other implications on the Galactic halo 
objects, damped Lya (DLA) systems and cosmology. 

2. TYPE lA SUPERNOVA PROGENITOR SYSTEM 

Our SD scenario has two progenitor systems: One 
is a red-giant (RG) companion with the initial mass of 
-^RG,o ~ ^Mq and the orbital period of tens to hun- 
dreds days (HKN96; HKN98). The other is a near 
main-sequence (MS) companion with the initial mass of 
-^MS,o ~ 2 — 3Mq and the period of several t enths of a 
day to several days (Li & van den Heuvel 1997; HKN98). 
In these SD scenarios, optically thick winds from the mass 
accreting WD play an essential role in stabilizing the mass 
transfer and escaping from forming a common envelope. 
The optically thick winds are driven b y a strong peak of 
OPA L opacity at logT(K) ~ 5.2 (e.g., Iglesias & Rogers 
1993). Since the peak is due to iron lines, the optically 



thick winds depend strongly on the metallicity (HK98). 

Figure 1 shows the metallicity dependence of the opti- 
cally thick winds. The strong winds are possible only for 
the region above the dashed line. The term "weak" im- 
plies that the wind velocity at the photosphere does not 
exceed the escape velocity there, that is, it cannot blow 
the accreted matter off the WD. For the metallicity as 



small 'AS Z = 0.001, the opacity peak at logT ^ 5.2 is 
very weak, being smaller than the peak of helium lines at 
logT ~ 4.6. Then, the wind is driven by the helium line 
peak rather than the iron line peak, which we call "He 
wind" instead of "Fe wind". Since only the ini tial WD 
mass of A fwn.n < 1.2M0 can produce an SN la (Nomotc 



fc Kondo 1991 ), SN la events occur only for the progeni- 
tors with [Fe/H] > —1.1, which is adopted in our chemical 
evolution model. 

Figure 2 shows the SN la regions in the diagram of the 
initial orbital period vs. the initial mass of the companion 
star for the initial WD mass of A/wd.o = 1-OMq (see HK98 
for other Mwd.o)- In these regions, the accretion from 
the companion star increases the WD mass successfully 
through the occurrence of SN la. The dashed line shows 
the case of solar abundance {Z = 0.02), while the solid line 
shows the much lower metallicity case of Z = 0.004. The 
size of these regions clearly demonstrate the metallicity ef- 
fect, i.e., SN la regions are much smaller for smaller metal- 
licity. The initial mass ranges of the companion stars for 
Z = 0.004 are O.9M0 < Mrg.o < 1.5Mq for the WD-f RG 
system and I.SMq < Mms,o < 2.6Mq for the WD-I-MS 
system. 



3. THE CHEMICAL EVOLUTION IN THE SOLAR 
NEIGHBORHOOD 

We use the chemical evolution model which allows the 
material infall from outside the disk region. For the in- 
fall rate, we adopt a formula which is proportional to 
texp[— t/r] with a infall timescale of r = 5 Gyr (Yoshii 
et al. 1996). The Galactic age is assumed to be 15 Gyr. 
The SFR is assumed to be proportional to the gas fraction 
with a constant rate coefficient of 0.37 Gyr~^. For the 
initial mass function (IMF), we assume a power-law mass 
spectrum with a Salpeter sl ope of a: = 1.35 in the range 
of O.O5M0 < M < 50Mq ([Tsujimoto et al. 1997]). We 



take the nucleosynthesis yields of SNe la and II from Tsu- 
jimoto et al. (1995; see Nomoto et al. 1997b, c for details) 
and the metallicity dependent main-sequence lifetime from 
Kodama (1997). 

For the SD scenario, the lifetime of SNe la is determined 
from the main-sequence lifetime of the companion star. 
We adopt the initial mass ranges of the binary companion 
stars obtained in section 2. The distribution function of 
the initial mass of the compa nions is taken from the mas s 



ratio distribution in binaries (Duguennoy & Mayor 1991 



which is approximated by a power-law mass spectrum with 
a slope X = 0.35. The fraction of primary stars of 3 — 8Mq 
which eventually produce SNe la is set to be 0.04 for both 
the WD-hMS and the WD-hRG systems, adjusted to re- 
produce the chemical evolution in the solar neighborhood. 
For the DD scenario we adopt the distribution function 
of the lifetime of SNe la by Tutukov & Yungelson (1994), 
majority of which is ~ 0.1 — 0.3 Gyr. 

Figure 3 shows the evolutionary change in [0/Fe] for 
three SN la models. The dotted line is for the DD sce- 
nario. The other lines are for our SD scenario with (solid 
line) and without (dashed line) the metallicity effect on 
SNe la. The results are as follows. 

• In the DD scenario the lifetime of the majority of 
SNe la is shorter than 0.3 Gyr. Then the decrease in 
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[0/Fe] starts at [Fe/H] ^ —2, which is too early com- 
pared with the observed decrease in [0/Fe] starting 
at [Fe/H] 1. 

• For the SD scenario with no metalhcity effect, the 
companion star with M ~ 2.6Mq evolves off the 
main-sequence to give rise to SNe la at the age of 
~ 0.6 Gyr. The resultant decrease in [0/Fe] starts 
too early to be compatible with the observations. 

• For the metallicity dependent SD scenario, SNe la 
occur at [Fe/H] > — 1, which naturally reproduce 

the observed break in [0/Fe] at [Fe/H] ~ — 1. 

We also perform the Kolmogorov-Smirnov (KS) test to 
estimate how well the p redicted abundance patt erns agree 
with the observed one ( Edvardsson et al. 1993| ), and tab- 
ulate the resultant probabilities in Table 1. The KS test 
is applied to the abundance distribution function in the 
range of —1.15 < [Fe/H] < 0.45 and to the evolutionary 
behavior of [0/Fe] in the range of -0.85 < [Fe/H] < 0.02. 
The predicted abundance distribution functions are identi- 
cal with the observational data with more than 84% prob- 
ability for both the DD and SD scenarios. For our SN 
la model with the metallicity effect, the calculated [0/Fe] 
against [Fe/H] fits well to the observations with 85% prob- 
ability. But when we neglect the metallicity effect, the 
probability decreases to 57%. For the DD scenario the 
calculated [0/Fe] is rejected with 82% probability. 



4. COSMIC SUPERNOVA RATE 

We apply our SN la progenitor scenario to predict 
the cosmic supernova rate history and the cosmic chem- 
i cal evolution corre s ponding to the obse rved cosmic SFR 
(iMadauetal. 1996| ; [Connolly ct al. 1997|) . The photomet- 
ric evolution is calculated with the spectral synthesis pop- 
ulation database taken from Kodama (1997). We adopt 
Ho = 50 km s"^ Mpc-\ = 0.2, Aq = 0, and the red- 
shift at the formation epoch of galaxies Z[ — 5. We use 
the initial c omoving density of gas ilgoo — 2 x 10^^ ( |Pe]| 
fc FaU 1995| ). 



Figure 4 shows the cosmic supernova rate per IO^^Lq 
per century (SNu). The long dashed line is for SNe II and 
the other lines for SNe la with (solid line) and without 
(dashed line) the metallicity effect. If we do not include 
the metallicity effect, the SN la rate is almost flat from 
the present to higher redshift, and decreases toward the 
formation epoch of galaxies. If we include the metallic- 
ity effect, the SN la rate drops at z ~ 1.2, where the iron 
abundance of the gas in the universe is too low (i.e., [Fe/H] 
< — 1) for the progenitors of SNe la to make explosions. 
The redshift where the SN la rate drops is determined by 
the speed of the chemical enrichment, which d epends on 
the effec t of dust extinction on the cosmic SFR (Pettini et 
al. 1997), cosmology, galaxy formation epoch and the ini- 
tial gas density. Taking into account these uncertainties, 
the break in the SN la rate occurs at z = 1 — 2. We should 
note that galaxies being responsible for the global SFR 
have different heavy-element enrichment timescale, some 
of which are starburst galaxies and ellipticals achieving 
[Fe/H] ~ —1 in much shorter timescale than in Figure 4. 
In such galaxies SNe la can occur even at z > 1 — 2. 



5. CONCLUSIONS AND DISCUSSION 

We introduce a metallicity dependence of the SN la rate 
in the Galactic and cosmic chemical evolution models. In 
our scenario involving a strong wind from WDs, little SNe 
la occur at [Fe/H] < — 1. Our model successfully repro- 
duces the observed chemical evolution in the solar neigh- 
borhood. We make the following predictions which can 
test this metallicity effect. 1) SNe la are not found in the 
low-iron abundance environments such as dwarf galaxies 
and the outskirts of spirals. 2) The cosmic SN la rate 
drops at z ~ 1 — 2 due to the low-iron abundance, which 
can be observed with the Next Generation Space Tele- 
scope. At z > 1 — 2, SNe la can be found only in the 
environments where the timescale of metal enrichment is 
sufficiently short as in starburst galaxies and ellipticals. 

The low-metallicity inhibition of SNe la can shed new 
light on the following issues: 

1) It imposes a limit on the metallicity range of the SN 
la progenitors. This would imply that "evolution effects" 
are relatively small for the use of high redshift SNe la to 
determine the cosmological para meters. 

2) Microlensing experiments (Alcock et al. 1997) sug- 



gest the WD-dominated Galactic halo. The existence of 
so many WDs results in too much iron enric hment from 
SNe la (Canal, Isern fc Ruiz-Lapuente 1997). However, 



WDs of halo populations are poor producers of SNe la, 
so that the iron abundance of the halo may not put the 
strong constraint on the contribution of WDs. Note, how- 
ever, that there exist seve ral other arguments against the 
WD-domina 1 
[Mould 1997|) 



1 a gaii 

WD-domina ted halo (e.g., [Chariot fc Silk 1995[ ; pibson & 



3) The Galactic halo is a low-metallicity system with 
[Fe/H] < — 1 and have an abundance pattern of gen- 
uine SN II origin, i.e., the overa bundances of g-elcmcnts 
relative to Fe as [a/Fe] > (e.g.. Wheeler, Sneden fc Tru-| 
ran 1989). However ther e exist age differences of several 
Gyrs amo ng the clusters ([Chaboyer, Demarque fc Saraje 



dini 1996) as well as field stars ( [Schuster fc Nissen 1989 } 



Since the shortest lifetime of SNe la is ~ 0.6 Gyr for the 
MS-f WD close binary systems, SN la contamination would 
be seen in [a/Fe] if there were no metallicity effect on SNe 
la. This apparent discrepancy between the age difference 
and the high [a/Fe] can be resolved by the low-metallicity 
inhibition of SNe la. 

4) Similar interpretation holds also for DLA systems. 
The DLA systems observed at 0.7 < z < 4.4 have [Fe/ H] = 
-2.5 to -1 and indicate [a/Fe] > ( |Lu ct al. 19961 ). Lu 
et al. (1996) suggested there may exist the agc-metallicity 
relation in DLA systems, which implies that DLA systems 
grown through a common chemical history spanning over 
several Gyrs. If so, the SN Il-like abundance pattern in 
DLA systems needs the introduction of the metallicity de- 
pendent SN la rate to avoid the contamination of SN la 
products. 
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Table 1 

Probabilities given by the Kolmogorov-Smirnov test 
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85.0% 89.5% 
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57.2% 84.8% 
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17.6% 86.9% 
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Fig. 1. — Mctallicity dependence of optically thick winds is shown in WD mass vs. mctallicity diagram. We regard the 
wind as "strong" if the wind velocity at the photosphere exceeds the escape velocity but "weak"' if the wind velocity is 
lower than the escape velocity. The term of "He" or "Fe" wind denotes that the wind is acxx^lerated by the peak of iron 
lines near logr(K) ~ 5.2 or of helium lines near logT(K) ~ 4.6. The dashed line indicates the demarcation between the 
"strong" wind and the "weak" wind. 
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Fig. 2. — The regions of SNe la is plotted in the initial orbital period vs. the initial companion mass diagram for the 
initial WD mass of Mwd,o = I-OMq. The dashed and solid lines represent the cases of solar abundance {Z = 0.02) and 
much lower mctallicity ot Z = 0.004, respectively. The left and the right regions correspond to the WD+MS and the 
WD+RG systems, respectively. 
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Fig. 3.— The evolutionary change in [0/Fe] against [Fe/H] for three SN la models. The dotted line is for the DD 
scenario where SNe la occur by merging of two WDs at a rate given by Tutukov & Yimgclson (1994). The other lines 
are for our SD scenario with (solid line) and without (dashed line) the mctallicity effect on SNe la. Observational data 
sources: filled circles, Edvardsson et al. (1993); open triangles, Barbuy & Erdelyi-Mendes (1989); stars, Nissen et al. 
(1994); open squares, Gratton (1991). 




Fig. 4. — The cosmic supernova rate per lO^^L© per century (SNu). The long-dashed line is for SNe II and the other 
lines for SNe la with (solid line) and without (dashed line) the metallicity effect. The filled circle is the observed SN la 
rate at .2 ~ 0.4 (Pain et al. 1996). The iron-abundance scale in the abscissa is calculated with the metallicity effect on 
SNe la. 



